SimArray: a user-friendly and user-configurable microarray design tool.

BACKGROUND: Microarrays were first developed to assess gene expression but are now also used to map protein-binding sites and to assess allelic variation between individuals. Regardless of the intended application, efficient production and appropriate array design are key determinants of experimental success. Inefficient production can make larger-scale studies prohibitively expensive, whereas poor array design makes normalisation and data analysis problematic. RESULTS: We have developed a user-friendly tool, SimArray, which generates a randomised spot layout, computes a maximum meta-grid area, and estimates the print time, in response to user-specified design decisions. Selected parameters include: the number of probes to be printed; the microtitre plate format; the printing pin configuration, and the achievable spot density. SimArray is compatible with all current robotic spotters that employ 96-, 384- or 1536-well microtitre plates, and can be configured to reflect most production environments. Print time and maximum meta-grid area estimates facilitate evaluation of each array design for its suitability. Randomisation of the spot layout facilitates correction of systematic biases by normalisation. CONCLUSION: SimArray is intended to help both established researchers and those new to the microarray field to develop microarray designs with randomised spot layouts that are compatible with their specific production environment. SimArray is an open-source program and is available from http://www.flychip.org.uk/SimArray/.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Adaptive evolution of Toll-like receptor 5 in domesticated mammals

<p>Abstract</p> <p>Background</p> <p>Previous studies have proposed that mammalian toll like receptors (TLRs) have evolved under diversifying selection due to their role in pathogen detection. To determine if this is the case, we examined the extent of adaptive evolution in the TLR5 gene in both individual species and defined clades of the mammalia.</p> <p>Results</p> <p>In support of previous studies, we find evidence of adaptive evolution of mammalian TLR5. However, we also show that TLR5 genes of domestic livestock have a concentration of single nucleotide polymorphisms suggesting a specific signature of adaptation. Using codon models of evolution we have identified a concentration of rapidly evolving codons within the TLR5 extracellular domain a site of interaction between host and the bacterial surface protein flagellin.</p> <p>Conclusions</p> <p>The results suggest that interactions between pathogen and host may be driving adaptive change in TLR5 by competition between species. In support of this, we have identified single nucleotide polymorphisms (SNP) in sheep and cattle TLR5 genes that are co-localised and co-incident with the predicted adaptive codons suggesting that adaptation in this region of the TLR5 gene is on-going in domestic species.</p

Saliva samples are a viable alternative to blood samples as a source of DNA for high throughput genotyping.

BACKGROUND: The increasing trend for incorporation of biological sample collection within clinical trials requires sample collection procedures which are convenient and acceptable for both patients and clinicians. This study investigated the feasibility of using saliva-extracted DNA in comparison to blood-derived DNA, across two genotyping platforms: Applied Biosystems Taqman™ and Illumina Beadchip™ genome-wide arrays. METHOD: Patients were recruited from the Pharmacogenetics of Breast Cancer Chemotherapy (PGSNPS) study. Paired blood and saliva samples were collected from 79 study participants. The Oragene DNA Self-Collection kit (DNAgenotek®) was used to collect and extract DNA from saliva. DNA from EDTA blood samples (median volume 8 ml) was extracted by Gen-Probe, Livingstone, UK. DNA yields, standard measures of DNA quality, genotype call rates and genotype concordance between paired, duplicated samples were assessed. RESULTS: Total DNA yields were lower from saliva (mean 24 μg, range 0.2-52 μg) than from blood (mean 210 μg, range 58-577 μg) and a 2-fold difference remained after adjusting for the volume of biological material collected. Protein contamination and DNA fragmentation measures were greater in saliva DNA. 78/79 saliva samples yielded sufficient DNA for use on Illumina Beadchip arrays and using Taqman assays. Four samples were randomly selected for genotyping in duplicate on the Illumina Beadchip arrays. All samples were genotyped using Taqman assays. DNA quality, as assessed by genotype call rates and genotype concordance between matched pairs of DNA was high (>97%) for each measure in both blood and saliva-derived DNA. CONCLUSION: We conclude that DNA from saliva and blood samples is comparable when genotyping using either Taqman assays or genome-wide chip arrays. Saliva sampling has the potential to increase participant recruitment within clinical trials, as well as reducing the resources and organisation required for multicentre sample collection.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Copynumber: Efficient algorithms for single- and multi-track copy number segmentation.

BACKGROUND: Cancer progression is associated with genomic instability and an accumulation of gains and losses of DNA. The growing variety of tools for measuring genomic copy numbers, including various types of array-CGH, SNP arrays and high-throughput sequencing, calls for a coherent framework offering unified and consistent handling of single- and multi-track segmentation problems. In addition, there is a demand for highly computationally efficient segmentation algorithms, due to the emergence of very high density scans of copy number. RESULTS: A comprehensive Bioconductor package for copy number analysis is presented. The package offers a unified framework for single sample, multi-sample and multi-track segmentation and is based on statistically sound penalized least squares principles. Conditional on the number of breakpoints, the estimates are optimal in the least squares sense. A novel and computationally highly efficient algorithm is proposed that utilizes vector-based operations in R. Three case studies are presented. CONCLUSIONS: The R package copynumber is a software suite for segmentation of single- and multi-track copy number data using algorithms based on coherent least squares principles.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Genomic and protein expression analysis reveals flap endonuclease 1 (FEN1) as a key biomarker in breast and ovarian cancer

FEN1 has key roles in Okazaki fragment maturation during replication, long patch base excision repair, rescue of stalled replication forks, maintenance of telomere stability and apoptosis. FEN1 may be dysregulated in breast and ovarian cancers and have clinicopathological significance in patients. We comprehensively investigated FEN1 mRNA expression in multiple cohorts of breast cancer [training set (128), test set (249), external validation (1952)]. FEN1 protein expression was evaluated in 568 oestrogen receptor (ER) negative breast cancers, 894 ER positive breast cancers and 156 ovarian epithelial cancers. FEN1 mRNA overexpression was highly significantly associated with high grade (p= 4.89 x 10 - 57) , high mitotic index (p= 5.25 x 10 - 28), pleomorphism (p= 6.31 x 10-19), ER negative (p= 9.02 x 10-35 ), PR negative (p= 9.24 x 10-24 ), triple negative phenotype (p= 6.67 x 10-21) , PAM50.Her2 (p=5.19 x 10-13 ), PAM50.Basal (p=2.7 x 10-41), PAM50.LumB (p=1.56 x 10-26), integrative molecular cluster 1 (intClust.1) ( p=7.47 x 10-12), intClust.5 (p=4.05 x 10-12) and intClust. 10 (p=7.59 x 10-38 ) breast cancers. FEN1 mRNA overexpression is associated with poor breast cancer specific survival in univariate (p=4.4 x 10-16) and multivariate analysis (p=9.19 x 10-7). At the protein level, in ER positive tumours , FEN1 overexpression remains significantly linked to high grade, high mitotic index and pleomorphism (ps< 0.01). In ER negative tumours, high FEN1 is significantly associated with pleomorphism, tumour type, lymphovascular invasion, triple negative phenotype, EGFR and HER2 expression (ps<0.05). In ER positive as well as in ER negative tumours, FEN1 protein over expression is associated with poor survival in univariate and multivariate analysis (ps<0.01). In ovarian epithelial cancers , similarly, FEN1 overexpression is associated with high grade, high stage and poor survival (ps<0.05). We conclude that FEN1 is a promising biomarker in breast and ovarian epithelial cancer

An open-access long oligonucleotide microarray resource for analysis of the human and mouse transcriptomes

Two collections of oligonucleotides have been designed for preparing pangenomic human and mouse microarrays. A total of 148 993 and 121 703 oligonucleotides were designed against human and mouse transcripts. Quality scores were created in order to select 25 342 human and 24 109 mouse oligonucleotides. They correspond to: (i) a BLAST-specificity score; (ii) the number of expressed sequence tags matching each probe; (iii) the distance to the 3′ end of the target mRNA. Scores were also used to compare in silico the two microarrays with commercial microarrays. The sets described here, called RNG/MRC collections, appear at least as specific and sensitive as those from the commercial platforms. The RNG/MRC collections have now been used by an Anglo-French consortium to distribute more than 3500 microarrays to the academic community. Ad hoc identification of tissue-specific transcripts and a ∼80% correlation with hybridizations performed on Affymetrix GeneChip™ suggest that the RNG/MRC microarrays perform well. This work provides a comprehensive open resource for investigators working on human and mouse transcriptomes, as well as a generic method to generate new microarray collections in other organisms. All information related to these probes, as well as additional information about commercial microarrays have been stored in a freely-accessible database called MEDIANTE

The Early Effects of Rapid Androgen Deprivation on Human Prostate Cancer.

The androgen receptor (AR) is the dominant growth factor in prostate cancer (PCa). Therefore, understanding how ARs regulate the human transcriptome is of paramount importance. The early effects of castration on human PCa have not previously been studied 27 patients medically castrated with degarelix 7 d before radical prostatectomy. We used mass spectrometry, immunohistochemistry, and gene expression array (validated by reverse transcription-polymerase chain reaction) to compare resected tumour with matched, controlled, untreated PCa tissue. All patients had levels of serum androgen, with reduced levels of intraprostatic androgen at prostatectomy. We observed differential expression of known androgen-regulated genes (TMPRSS2, KLK3, CAMKK2, FKBP5). We identified 749 genes downregulated and 908 genes upregulated following castration. AR regulation of α-methylacyl-CoA racemase expression and three other genes (FAM129A, RAB27A, and KIAA0101) was confirmed. Upregulation of oestrogen receptor 1 (ESR1) expression was observed in malignant epithelia and was associated with differential expression of ESR1-regulated genes and correlated with proliferation (Ki-67 expression).We thank CRUK, The NIHR, The Academy of Medical Sciences(RG:63397) and the National Cancer Research Prostate Cancer: Mechanisms of Progression and Treatment (ProMPT) collaborative (G0500966/75466), Hutchison Whampoa Limited, the Human Research Tissue Bank (Addenbrooke’s Hospital, supported by the NIHR Cambridge BRC), and Cancer Research UK

HNF1B variants associate with promoter methylation and regulate gene networks activated in prostate and ovarian cancer

Two independent regions within HNF1B are consistently identified in prostate and ovarian cancer genome-wide association studies (GWAS); their functional roles are unclear. We link prostate cancer (PC) risk SNPs rs11649743 and rs3760511 with elevated HNF1B gene expression and allele-specific epigenetic silencing, and outline a mechanism by which common risk variants could effect functional changes that increase disease risk: functional assays suggest that HNF1B is a pro-differentiation factor that suppresses epithelial-to-mesenchymal transition (EMT) in unmethylated, healthy tissues. This tumor-suppressor activity is lost when HNF1B is silenced by promoter methylation in the progression to PC. Epigenetic inactivation of HNF1B in ovarian cancer also associates with known risk SNPs, with a similar impact on EMT. This represents one of the first comprehensive studies into the pleiotropic role of a GWAS-associated transcription factor across distinct cancer types, and is the first to describe a conserved role for a multi-cancer genetic risk factor

SPAG5 as a prognostic biomarker and chemotherapy sensitivity predictor in breast cancer: a retrospective, integrated genomic, transcriptomic, and protein analysis

Background: Although the use of proliferation markers/profiles has been recommended when choosing the appropriate systemic-treatment for breast cancer (BC), the best molecular-marker/test that should be used needs to be identified. Methods: To identify factors that drive proliferation and its associated features in BC an artificial neural network (ANN) based integrative data-mining methodology was applied to three cohorts [(Nottingham-discovery (ND), Uppsala and METABRIC (Molecular Taxonomy of Breast Cancer International Consortium)]. The most prominent genes in the resulting interactome-map were then identified. Given that SPAG5 was associated with many features of proliferation, featured prominently in the interactome-map and has a fundamental role in mitotic-progression,, we hypothesized that it could be a better indicator of proliferation activity. (BC). Subsequently to test if it could provide a more accurate guide for the delivery of systemic therapies in BC, we investigated the clinico-pathological utility of SPAG5: gene copy number aberrations (CNAs); mRNA and protein expression, in over 10,000 BCs. Integrated analysis of SPAG5-gene CNAs, transcript and protein expression was conducted in the ND cohort (n=171) and validated in the METABRIC cohort (n=1980). In addition, the associations of SPAG5 CNAs, transcript and/or protein with BC specific survival (BCSS), disease free survival (DFS) and/or distant relapse free survival (DRFS) were analysed in multiple cohorts including Uppsala (n=249), METABRIC, three-untreated lymph node (LN) negative cohorts (n=684), a combined multicentre clinical data set (n=5439), Nottingham historical early-stage-primary BC (Nottingham-HES-BC; n=1650), Nottingham oestrogen receptor (ER) negative BC (n=697), Nottingham anthracycline-Neoadjuvant-chemotherapy (Nottingham-AC-Neo-ACT; n=200), and MD Anderson Cancer Centre Taxane/anthracycline (MDACC-T/AC-Neo-ACT; n=508) cohorts. The association of SPAG5 transcript and protein expression with pathological response rate (pCR) were also tested in [MDACC-T/AC-Neo-ACT (n=508) and the phase II trial NCT00455533; n=253)] and [Nottingham-AC-Neo-ACT (n=200)] cohorts; respectively. Findings: SPAG5 gene gain/amplification at the Ch17q11·2 locus was found in 10.4% of BC (206/1980 (; METABRIC) and was reported in 19·4% of PAM50-HER2 (46/237) and 17·8% of PAM50-LumB (87/488). SPAG5-CNA gain/amplification and high SPAG5-transcript and SPAG5-protein were associated with increased risk of death from BC [Uppsala; (HR (CI 95%): 1·50 (1·18-1·92); p=0·00010, METABRIC; (HR (CI 95%): 1·68 (1·40-2·01) p<0·0001), and Nottingham-HSE-BC; (HR (CI 95%): 1·68 (1·32-2·12), p<0·0001); respectively]. Multivariable Cox regression models, including other validated-prognostic factors, (Uppsala: age, size, LN-stage, genomic grade index (GGI), ER, TP53 mutation and MKi67; METABRIC: age, size, LN-stage, histologic-grade, ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), hormone-therapy, chemotherapy, interaction term of SPAG5 and both chemotherapyy and hormonotherapy; Desmedt-untreated LN- cohort: ER, Nottingham prognostic index (NPI), 76-gene prognostic signature (Veridex) and Adjuvant-Online (AOL); Nottingham-HES-BC: menopausal status, size, LN- stage, histologic-grade, ER, PR, HER2, ki67, hormone-therapy, chemotherapy, interaction term of SPAG5 and both chemotherapy[y and hormonotherapy), showed that high SPAG5-transcript and high SPAG5-protein were associated with shorter BCSS [Uppsala: (HR (CI 95%): 1·62 (1·03-2·53) p=0·036); METABRIC: (HR (CI 95%): 1·27 (1·02-1·58) p=0·034); Desmedt-untreated LN- cohort: (HR (CI 95%): 2·34 (1·24-4·42) p=0·0090), and Nottingham-HES-BC (HR (CI 95%): 1·73 (1·23-2·46) p=0·0020); respectively]. In ER-negative-BC with high SPAG5-protein, administration of anthracycline-adjuvant-chemotherapy had reduced the risk of death by 60% compared to chemotherapy-naive (HR (95% CI): 0·37 (0·20-0·60); p=0·0010). A multivariable Cox regression analysis, which included other validated prognostic factors for chemotherapy (e.g., menopausal status, size, lymph node stage, histologic grade, ER, PR, HER2, Bcl2, chemotherapy, interaction term of SPAG5 and both chemotherapy[y), revealed that SPAG5-transcript+ was independently associated with decreased risk of DRFS after receiving Taxane/anthracycline-Neo-ACT [MDACC-T/AC-Neo-ACT: (HR (CI 95%): 0·68 (0·48-0·97); p=0·0070)]. In multivariable logistic regression analysis, both SPAG5-transcript+ and SPAG5-protein+ and were independent predictors for higher pCR after combination-cytotoxic chemotherapy [MDACC-T/AC-Neo-ACT: (OR (95% CI) 1·71 (1·07-2·74); p=0·024), and Nottingham-AC-Neo-AC: (OR (95% CI): 8·75 (2·42-31); p=0·0010); respectively]. Interpretation: SPAG5 is a novel amplified gene on Ch17q11.2 in PAM50-LumB and PAM-HER2 BC, and its transcript and protein products are independent prognostic and predictive biomarkers, with potential clinical utility as a biomarker for combination cytotoxic chemotherapy sensitivity, especially in ER- BC

The somatic mutation profiles of 2,433 breast cancers refines their genomic and transcriptomic landscapes

The genomic landscape of breast cancer is complex, and inter- and intra-tumour heterogeneity are important challenges in treating the disease. In this study, we sequence 173 genes in 2,433 primary breast tumours that have copy number aberration (CNA), gene expression and long-term clinical follow-up data. We identify 40 mutation-driver (Mut-driver) genes, and determine associations between mutations, driver CNA profiles, clinical-pathological parameters and survival. We assess the clonal states of Mut-driver mutations, and estimate levels of intra-tumour heterogeneity using mutant-allele fractions. Associations between PIK3CA mutations and reduced survival are identified in three subgroups of ER-positive cancer (defined by amplification of 17q23, 11q13-14 or 8q24). High levels of intra-tumour heterogeneity are in general associated with a worse outcome, but highly aggressive tumours with 11q13-14 amplification have low levels of intra-tumour heterogeneity. These results emphasize the importance of genome-based stratification of breast cancer, and have important implications for designing therapeutic strategies.The METABRIC project was funded by Cancer Research UK, the British Columbia Cancer Foundation and Canadian Breast Cancer Foundation BC/Yukon. This sequencing project was funded by CRUK grant C507/A16278 and Illumina UK performed all the sequencing. The authors also acknowledge the support of the University of Cambridge, Hutchinson Whampoa, the NIHR Cambridge Biomedical Research Centre, the Cambridge Experimental Cancer Medicine Centre, the Centre for Translational Genomics (CTAG) Vancouver and the BCCA Breast Cancer Outcomes Unit. We thank the Genomics, Histopathology, and Biorepository Core Facilities at the Cancer Research UK Cambridge Institute, and the Addenbrooke’s Human Research Tissue Bank (supported by the National Institute for Health Research Cambridge Biomedical Research Centre).This is the final version of the article. It first appeared from Nature Publishing Group via http://dx.doi.org/10.1038/ncomms1147